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INTRODUCTION 


Robot dynamics and control will become an important issue for implementing productive 
platforms in Space. Robotic operations will become necessary for man-tended stations and for 
efficient performance of routine operations in a manned platform. The current constraints on the 
use of robotic devices in a microgravity environment appears to be due to an anticipated increase 
in acceleration levels due to manipulator motion and for safety concerns. The objective of this 
study will be to provide baseline data to meet that need. 

Most texts and papers dealing with kinematics and dynamics of robots assume that the 
manipulator is composed of joints separated by rigid links. However, in recent years several 
groups have begun to study the dynamics of flexible manipulators, primarily for applying robots 
in space and for improving the efficiency and precision of robotic systems. Some of the reasons 
for pursuing research in this area are: 

1. Current industrial robot systems weigh substantially more than the payload that they 
are able to transport. For instance a 250 pound robot can usually only transport 8- 10 pounds at 
the end-effector. Hence by increasing the ratio of payload to arm weight will substantially 
improve the economics for using robots through more efficient usage of materials and power. 

2. Improving the precision of robotic systems for performing delicate tasks will enable a 
broader base for robotics applications. 

3. A number of applications, such as welding, deburring, etc. require trajectory tracking 
which may exceed the capability of today's industrial robot base, particularly when unexpected 
parameters cause the robot motion to deviate from the pre-programmed paths. 

SPACE ROBOTICS 

Robotic systems which are being planned for implementation in space have a number of 
constraints to overcome. 1 ‘ 4 The need to quantify some of the characteristics that are required 
for structures such as the Space Station is very important at this time, even though the IOC will 
be 5 -10 years from now. 
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Additional concepts which have to be worked out in any robotic implementation for a 
space platform include teleoperation and degree of autonomous control. Lightweight arms are 
necessary for space, primarily for the cost benefits derived from their reduced weight. However, 
lighter weight arms have to necessarily flex during movement. Flexure of the arms performing a 
task requiring precision requires some control mechanism to insure that the end-effector is at the 
proper place and orientation with respect to the workpiece during any contact period. 

Flexing motions of the arm can cause (1) accelerations to feed back into the base support 
of the robot and into the Space Station structure, (2) transmit accelerations into the sample being 
transported, or (3)take forever to perform a task. The first effect is obviously detrimental to the 
microgravity environment of the Space Station, while the second will impact experiments such 
as the delicate protein crystals which are to be grown in space and probably transported with 
robotic or telerobotic arms. In some cases slow movements may be acceptable for (3); however 
it certainly will not be suitable for most tasks required for operation on the Space Station. 

One must also include the reasoning that for man and robots to co-exist in the space 
environment, the robot must be non-threatening to man. Lightweight arms satisfy that criteria. 
For space applications, a Cincinnati-Milacron T3 is not only over-weight, but it also is 
threatening to humans trespassing in its working volume. 

PREVIOUS WORK 

The robot used for the initial study in NAS8-36620 was a UMI RTX robot, which was 
adapted to operate in a materials processing work cell to simulate sample changing in a 
microgravity environment. The robotic workcell was flown several times on the KC-135 aircraft 
at Ellington Field. The primary objective of the initial flights was to determine operating 
characteristics of both the robot and the operator in the variable gravity of the KC-135 during 
parabolic manuevers. 

This study demonstrated that the KC-135 aircraft can be used for observing dynamics of 
robotic manipulators. We also observed the difficulties associated with humans performing 
teleoperation tasks during varying G levels and provided insight into some areas in which the use 
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of artificial techniques would provide improved system performance. 

Allowing some degree of autonomy due to time delay communications for teleoperation 
over large distances is also necessary. The robot controls then will have certain motions 
embedded in the control software that do not need explicit operator communication, except for 
abort. Time delays make for precarious situations in performing teleoperation from large 
distances. Tasks can then be accomplished in a more reasonable manner and more successfully. 

A major goal of the earlier study was to evaluate a small robot system, such as the UMI 
RTX, for materials processing applications in low gravity and to determine the characteristics of 
a robot arm in a space environment, particularly with respect to accelerations which might 
impact materials grown on a space platform. 

A materials transfer workcell was assembled to simulate the changing of sample 
ampoules as might occur aboard a space laboratory. Accelerometer packages were included for 
determining the G levels of the workcell and the at the end-effector. Several KC-135 flights 
were made with the workcell, improving some data taking capabilities each mission. The human 
operator was able to train the robot to perform a materials transfer function within the 20 
seconds desired. The first computer used with the experiment was a Toshiba 3100 which did not 
allow for both control of the robot and reading of the accelerometer pacKage at the same time. 
The microprocessor in that system was an 80286. In order to try to improve upon the multi- 
tasking capability of the computer, a 386 based computer was selected. This choice enabled 
some improvements in the data acquisition process; but the multi-tasking software used at that 
time still did not permit the I/O commands to the robot to operate properly. Consequently we 
never did get to control the robot and take acceleration data simultaneously in these flights. 

A number of lessons were learned with this series of experiments. The RTX robot uses 
plastic belts for actuation of the links and optical encoders for position and velocity control. The 
slippage and flexing of the belts caused excessive jitter and accelerations at the end-effector. We 
believe that the belt-driven actuation would not be acceptable for experiments such as the protein 
crystal growth studies due to the lack of control of accelerations at the end-effector. The control 
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system; however, is PID and appeared to work well whether the task was learned in 1 G and 
performed in low G or vice-versa. However, it was tedious to teach the robot during parabolas, 
mainly because we had few visual aids to assist in the correct orientation of the end effector. A 
particularly sensitive task was inserting the sample ampoule into its holder. A borrowed fiber- 
optic borescope provided little depth perception and was not useful for this study. In addition 
teaching a robotic device for precision movements can certainly be improved through more 
innovative approaches using embedded sensors or vision systems with some autonomous local 
control using artificial intelligence techniques. 

For example, protein crystal growth experimentation has great interest at this time and 
will be a primary space experiment. The robot which has historically been used within that 
program for sample preparation and is within a reasonable price range, is the Zymate II robot 
manufactured by the Zymark Corporation. It possesses a more sophisticated controller than the 
UMI RTX used in the earlier study and has a greater potential for working in a multi-tasking 
environment. 

The computer controller is important because robot programming, robot execution and 
data acquisition of the accelerometers has to be performed during parabolic manuevers. The MS- 
DOS platform used in earlier study doesn't allow all those functions to occur simultaneously. We 
chose to move to a more natural multi-tasking environment offered by the Commodore Amiga 
for the second phase activities. 

The focus on protein crystal growth also required a different robot workcell to be 
designed and fabricated. Another major concept which might be important in terms of promoting 
telescience experiments is to use the KC-135 to implement the above experiments with remote 
manipulation from the ground. The KC-135 aircraft facility personnel at JSC have indicated 
interest in these types of experiments using the TDDRS satellite for transmission purposes. 

CONCEPTS OF ROBOT DYNAMICS 

The Newton-Euler Equations state be used to describe the effective forces distributed 
throughout the robotic system. In this work, we are particularly interested in the residual forces 
present at the end-effector. 
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f i =1 m i r i 


Forces due to linear accelerations 


Tj - £ (IjWj + wj x IjWj) Torques due to rotational velocities and 

accelerations 

Note that these torques are referred to as dynamic torques since they arise only during motion, 
whereas static torques can exist such as when power is required to maintain a load against 
gravity. The three types of dynamic torques present in a moving robot arm are: 

inertial torques - proportional to joint accelerations, 
centripetal torques - proportional to velocities squared 
Coriolis torques - proportional to joint velocities from two 
different links. 

The inertial forces and torques arise from the normal action/reaction forces when 
accelerating a body. Obviously in order to move the arm from one resting position to another, 
both acceleration and decceleration will occur. Centripetal torques arise from contrained rotation 
about an axis. These torques will be present whenever a single rotary motion is executed. The 
Coriolis forces arise from the interaction of two simultaneously rotating systems. 

When expanded in this manner, the dynamic equations increase in complexity with the 
number of joints. Each manipulator will have its own characteristic set of relations describing 
these torques. The complexity of the relations vary due to the joint interaction of these inertial, 
centripetal, and Coriolis torques. 

In order to provide some control over a robotic system in terms of a consistent and 
reproducible path generation for varying payloads, a PID controller is normally used. Figure 1 
shows in a simplified form the normal block diagram used for such a controller.^ In order to 
compensate for gravity, for example, the PID controller allows torque to be delivered to the 
motors just to maintain position. Additional torque is required for movement. 
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A more complex controller is available in Treetops which includes more variables such 
as friction, gear backlash, and stiffness. The block diagram for Motor4 is shown in Figure 2. 

Figure 2. Block diagram for Motor4 as taken from the Treetops Users Manual 









RESULTS 


The work performed in this contract consisted of two distinct activities. One set was 
primarily concerned with the construction of a Zymate robotic workcell which could fly on the 
KC- 1 35 and take accelerometer data for determining the g-levels which would be present within 
the robotic system when it was operating in a reduced gravity environment. The other set of 
activities concentrated on ground-based applications in order to determine what improvements 
were able to be attaiined in improving the preciseness of performing delicate robotic tasks, even 
in 1-g. These were primarily performed at the Alabama Center for Advanced Technology 
Transfer, (ACATT) and were performed with a Puma 560 robot owned by the Boeing Aerospace 
Corporation. 

ZYMATE ROBOT CELL 

The robot workcell which was constructed for the KC-135 simulations is shown in 
Figure 3. The workcell was flown in June, 1990 for a systems check-out flight. Unfortunately 
during the remainder of the contract, the KC-135 did not fly so we were unable to finish up all 
the sub-systems performing in the reduced gravity environment. 

FIGURE 3. Schematic of Zymate Workcell for KC-135 Experiments 
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Figure 4. Acceration Measurements on Puma 560 Robot 

a.) Overlays of accelerations in x, y, and z vectors for a one rotational joint motion. 



Overlays of x, y, and z accelerations during multiple motion trajectories. 
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TREETOPS SIMULATIONS 

In parallel with the KC-135 experiments, we have also begun a number of simulated 
experiments with the intent of developing a better understanding of the dynamics of both rigid 
and flexible manipulators. The TREETOPS Simulation software was obtained from Dr. Henry 
Waites' group at Marshall Space Flight Center. One of the graduate students, Mr. Houchang Li 
has been performing these simulations on a VAX 785 located at the University. We feel that 
these simulations will be applicable to both reduced gravity and ground-based robotic 
applications. 

TREETOPS provides availability to a number of devices, actuators, sensors, etc. Figure 5 
shows a table listing some of more significant components used for robotics simulations. 


O * - 

COMPONENTS AVAILBLE TO THE USER FOR TREETOPS SIMULATION 

BODY 

HINGES 

SENSORS 

ACTUATORS 

FUNCTION 

GENERATORS 

INTERCONNECTIONS 





STEP 

LICARSPRNG 

RIGD 

ROTATION 

RATE GYRO 

REACTION JET 





QUADRATIC SPRNG 

FlEXBLE 

TRANSLATION 

RESOLVER 

HYDflAULC 

RAM 3 


MOON. 

ANGULAR 

CYLH3ER 

PULSE 

SOLD DAMPER 


NASTOAN 


MOMENT 





ACCELERATE* 

ACTUATOR 

TRIANGLE 

UPPER HARDSTOP 



VELOCITY 

TORQUE 

SfC 

LOWER HAROSTOP 



POSITION 

MOTOR BRAKE 

NOSE 




TACHOMETER 

MAGNET 

DOUBLET 




NTEGRATING 

GWBALCMG 

USER 



GYRO 

MU SENSOR 
POSfTDN 


VELOCITY 

VECTOR 

CMG RESOLVER 

CMG TACHOMETER 


REACT ON 
WHEEL 


MOTOR 

DRNES 
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In order to access these devices, TREETOPS provides several modules which allow one to enter 
the desired input parameters and simulate specific time histories. Figure 6 shows how these 
modules fit together. 


Figure 6. TREETOPS Architecture 




In cooperation with the Alabama Center for Advanced Technology Transfer, a graduate 
student in Mechanical Engineering has been using TREETOPS to model the dynamics of several 
industrial robot systems. Mr. Houchang Li has worked on models of both a Puma 560 for 
activities at ACATT to improve the the precision of robotic systems and the Zymate for this 
project. Examples of the dynamics which can be modeled using TREETOPS can be seen in the 
following analysis. 

Using the motor4 model available in TREETOPS to drive the robot arm in three different 
motions, the effect of gravity can be investigated. Three different trajectories are displayed in 
Figure 7: 

a. ) A vertical motion into the motor labelled 2. 

b. ) A rotation about the vertical axis by motor 1. 
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c.) A prismatic motion outward using motor 3. 

Figure 7. Illustration showing sequences of motion for Zymate robot used in TREETOPS 
simulation. 



The three trajectories are entered into TREETOPS as motions of the three links and the 
end-effector and the resulting time history of the motions are plotted in Figure 8. Note that the 
vertical motion follows the input pulse into motor 2 for the vertical motion from 0 to 2.5 
seconds. This trace is labelled 1 in the figure. In order to simplify the chart, the other two input 
waveforms are left off; however all three motions of the robot are plotted. 
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Figure 8. TREETOPS plot showing the motion of the end-effector for the trajectory described 
in Figure 7. 



The corresponding motions for the other links and the end-effector can be determined by 
relating the time motion that is occuring with the motor activated waveforms from Figure 7. 
Note that the label 2 refers to a displacement about the rotation axis, the label 3 refers to a 
displacement along the prismatic axis, and 4 is the vertical displacement of the end-effector. 

The effect of gravity can be seen in the next two figures. Figure 9 shows the how the 
dynamics (or accelerations) at the end-effector occur during the performance of the trajectories 
with the gravity vector turned on. In Figure 10, the same motions are performed with the gravity 
vector turned off. Note that in the absence of gravity, there are more acceleration components 
occuring at the end-effector than in 1 g. Gravity does then tend to damp out some oscillations of 
the end-effector. In microgravity then for the same manipulator and control scheme, one would 
expect more oscillations. 
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Figure 9. The TREETOPS output for the vertical direction for the simulation decribed above 


with gravity. 






Conclusions 

This work has resulted in a robot workcell which is qualified to fly on the KC-135 for 
measuring the responses of robotic devices in low gravity. Due to the 20-25 seconds of low 
gravity available during parabolic flight, one has to define trajectories which fit into that time 
period. In simulating sample changing operations, that type of trajectory works very well. Other 
trajectories should be suitable; however, if they are functional tasks that a robotic system would 
be performing in space, then they are more easily accepted as useful trajectories. More woerk 
needs to be done in identifying such tasks. 

Another improvement that would be beneficial in the future is a better accelerometer 
system and interconnections to the computer. For this effort, accelerometers were borrowed 
from another KC-135 experiment. A smaller, lightweight triaxial accelerometer mount needs to 
replace the current single axis one used in this work due its overall mass. 

Other improvements conceptually include a wireless signal for the accelerometer 
measurement, so that the cables do not interfer with the motion of the robot arm and also to 
reduce the potential for the cables to impart vibration to the accelerometer itself. These ideas can 
be looked at more closely in future projects. 

During the course of this research a number of reports which dealt with acceleration 
control of robotic systems were found. The goals of these researchers were quite similar to ours 
in trying to improve precision and payload capability for industrial applications. These papers 
are listed in the bibliography and should be consulted for more information. 

Summary 

This research has generated some significant results in developing a robotic workcell for 
performing robotics research on the KC-135 in preparation for space-based robotics applications 
in the future. In addition, we have shown that TREETOPS can be used to simulate the dynamics 
of robot manipulators for both space and ground-based applications. 
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APPENDIX - Programs developed for the Amiga 2500 to acquire and process accelerometer 
data from the Zymzte robot. 


/************************************************************************ 
*readlc.c * 

* by David Gilliam UAH * 

* * 

* Reads one channel of the accelerometer and saves the data after * 

* averaging every group of four values. * 

* * 


******************************************************************** 


****^ 


#include <exec/types.h> 
#include <stdio.h> 
#include <p40.h> 
#include <dos.h> 


main() 

{ 

short data[24000]; 
short finaldata[6000]; 
short channel, success; 
short gain, rate, i; 
long samples; 
float volts = 0.0; 
int fh, error; 
char fname[40],c; 
unsigned int count; 

if((success = p40_init())==0) { 

printf("Couldn't configure the Proto-40k\n"); 
exit(0); 

} 

printf("\nlnput data file name : "); 
i = scanf("%s",fname); 
printf("\n"); 
c = getchar(); 

samples = 4096; 

gain = 0; /* gain of 10 */ 

channel = 0; 

rate = 400; /* in Hz */ 
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*p40_modeselect = 0; 

*p40_mux = (channel « 2); /* Set A/D channel */ 

/* p40_setgain(gain); 
p40_setrate(rate); */ 

printf("\nPress enter to begin taking data\n”); 
c = getchar(); 

printf("Taking data, please wait...\n"); 
if(p40_atod(data,rate,samples,0,gain) == 0) 
printf("Couldn't acquire data\n"); 
printf("Data acquisition completeAn"); 
printf("Smoothing data\n\n"); 

for(i=0;i<samples;i++) 

finaldatafi] = (data[4*i] + data[4*i+l] + data[4*i+2] + data[4*i+3])/4; 

printf("\nDisplay data? (y/n)"); 

c = getchar(); 

printf("\n"); 

if (c == V) { 

printf("Count Channel %02d\n", channel); 
printf(" An"); 

for (i=0; i<samples/4; i++) { 

volts = ((float)finaldata[i] / 204.8) - 10.0; 

printf("%04x %04d %7.4f\n", finaldata[i], finaldata[i], volts); 


} 

} /* end of if */ 

printf("\nSaving data to disk\n"); 

fh = creat(fname, 0); 

count = write(fh, &finaldata, 12000); 

error = close(fh); 


} 

/ ********* * ******** * ** * * 4c * * * * 4c * * * 4c * * * * * * * * * * * * * 4c * * * 4c 4c 4c * * * * * * * * * * * * * * * * * * 

* showlc.c * 

* by David Gilliam UAH * 


* Displays a 'single channel' data file * 

******* ****** * ********* * *** * * * ** * * * ** * * * * * * * * * * % * % * * * * * * * * * * * * * * * * * * * * * */ 


#include <exec/types.h> 
#include <stdio.h> 
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#include <p40.h> 
#include <dos.h> 
#include <fcntl.h> 


main() 

{ 

short data[6000]; 
short channel; 
short gain, rate, i; 
long samples; 
float volts = 0.0; 
int fh, error; 
char *fname,c; 
unsigned int count; 

printf("\nlnput data file name : "); 
i = scanf("%s",fname); 
printf(’V); 

c = getchar(); 

samples = 600; 

gain=l; /* gain of 10 */ 

channel = 0; 

rate = 100; /* in Hz */ 


} 


printf(”\n Loading data from disk\n"); 
fh = open(fname, 0_RDONLY, 0); 
count = read(fh, &data, 1 2000); 
error = close(fh); 


for (i=0; i<samples; i++) { 

volts = ((float)((data[i] / 204.8) - 10)); 


printf("%04x %04d %7.4f\n", data[i], data[i], volts); 


if(i % 21 == 20) { 
c = getchar(); 
if ( C == -q’ | c == 'Q') 
i = samples; 

} 


} 


^*************** ******* jic**************,^.^^^^^^^^^^^^^*^^^ 

* convl.c * 
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* by David Gilliam UAH * 


#include <exec/types.h> 

#include <stdio.h> 

#include <p40.h> 

#include <dos.h> 

#include <fcntl.h> 

main() 

{ 

FILE *fp; 
short data[6000]; 
short channel; 
short gain, rate, i; 
long samples; 
float volts = 0.0; 
int fh, error; 
char fname[40],c; 
unsigned int count; 

printf("\nlnput data file name : "); 
i = scanf("%s",fname); 
printf("\n"); 
c = getcharf); 

samples = 1200; 

gain = l; /* gain of 10 */ 

channel = 0; 

rate = 100; /* in Hz */ 

printf('\nLoading data from disk\n"); 
fh = openffname, O.RDONLY, 0); 
count = read(fh, &data, 1 2000); 
error = close(fh); 

strcat(fname,".asc"); 
fp = fopenffname, "w"); 
if (fp == 0) { 

printf("unable to open destination fileW); 
exit(); 

} 

/* printf("Count Channel 0\n", channel); 
printf(" - — -\n"); */ 
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fprintf(fp,'\naccelerometer data\n\n"); 
for (i=0; i<samples; i++) { 

/* volts = ((float)(data & OxOfff - 2048)/ 20480.0); */ 

volts = ((float) ((data[i] / 204.8) - 10)); 

/* volts = (volts - 0.6) / 0.6; */ 

fprintf(fp,"%7.4f\n", volts); 

) 

fcloseallO; 

} 
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